Creates the option to maintain resolution in the ocean. Used to resolve
fjords and bathymetry when coupling to MPAS-Ocean.

This commit creates the option to maintain high resolution in all ocean
cells (same resolution as what is prescribed at outlet glacier termini).
This is necessary to resolve the fjord geometry that thermal forcing
parameterizations rely on. To reduce time when generating these higher
resolution meshes, this commit also introduces an option to do an
initial cull of the bedMachine dataset using the extent of the
contintental shelf before generating the GIS mesh.

This commit moves the calculation of geom_points and geom_edges
using a geojson of the gis contintental shelf extent
to compass/landice/tests/greenland/mesh.py, instead of
compass/landice/mesh.py. If this options is not being used, thn
geom_points and geom_edges are calculated within build_cell_width, as
per usual

Co-authored-by: Trevor Hillebrand <[email protected]>

Re-introduces code that sets a large cell width for cells that will be
culled by distance from ice edge. Only applied when not using a geojson
to pre-cull domain

Ensures the bedTopography constraint on mesh resolution is enforced
regardless 'max_res_in_ocn'. Instead, grid resolution in the ocean
is maximized following the bedTopgraphy constraint, overwritting it
where necessary.

Introduces new landice test case to generate a domain for the
Uummannaq/Disko Bay region of GIS.

Create region masks for each tidewater outlet glacier in domain

Addresses bug introduced by the inclusion of the parallel_executable
argument in interp_gridded2mali

Adds surface and basal mass balances fields and their respective
uncertainties to the uummmannaq/disko mesh_gen test case. used for the
dhdt albany optimization